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Abstract 

We present a comprehensive optical study of the narrow gap FeSb2 semiconductor. From the 
optical reflectivity, measured from the far infrared up to the ultraviolet spectral range, we extract 
the complete absorption spectrum, represented by the real part a\{oj) of the complex optical 
conductivity. With decreasing temperature below 80 K, we find a progressive depletion of (J\(uj) 
below Eg ~ 280 cm -1 , the semiconducting optical gap. The suppressed (Drude) spectral weight 
within the gap is transferred at energies oj > E g and also partially piles up over a continuum of 
excitations extending in the spectral range between zero and E g . Moreover, the interaction of one 
phonon mode with this continuum leads to an asymmetric phonon shape. Even though several 
analogies between FeSb2 and FeSi were claimed and a Kondo-insulator scenario was also invoked 
for both systems, our data on FeSbi differ in several aspects from those of FeSi. The relevance 
of our findings with respect to the Kondo insulator description will be addressed. 

PACS numbers: 78.20.-e 
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I. INTRODUCTION 



FeSb2 represents an interesting case of narrow gap semiconductor, where a band of 
itinerant electron states originates in the d xy orbitals of the t 2g multiplet. The magnetic 
properties 1 ' 2 of FeSb 2 strongly resemble those seen in FeSi (Refs. 3 and 4). The magnetic 
susceptibility of FeSb2 shows a diamagnetic to paramagnetic crossover around 100 K (Ref. 
1); however in contrast to FeSi, with a very small low temperature (T) impurity tail in 
the diamagnetic region. The electrical transport in FeSb2 along two axes is semiconducting, 
whereas the third axis exhibits a metal-semiconductor crossover at temperature T cr ~ 40—80 
K, depending from the current alignment 1 . Furthermore, the thermal expansion and heat 
capacity of FeSb 2 at ambient pressure agrees with a picture of a temperature induced spin 
state transition within the Fe t 2g multiplet 2 . 

Aeppli and Fisk 5 proposed that the underlying physics of FeSi share common features 
with a class of rare earth compounds, known as hybridization gap semiconductors or Kondo 
insulators. The analogies between FeSi and FeSb 2 , pointed out above, would suggest that 
also FeSb 2 might belong to the class of ci-electron based Kondo insulators. Within the Kondo 
insulator's scenario, one considers two narrow hybridized bands of width W in the density 
of states, separated by the gap E g . At T = the electrons populate a lower hybridized 
band, and with the increase of T the electrons start populating the higher band, resulting 
in a thermally activated Pauli susceptibility 5 . This approach, successfully applied to FeSi, 
seems to work out for FeSb 2 , as well 1 ' 2 . This is of particular interest because of possible 
relations between theoretical descriptions of d- and /-electron systems. 

FeSi was also intensively investigated from the point of view of its absorption 
spectrum 6 ' 7,8 ' 9 . The charge excitation spectrum of FeSi is characterized at low temperatures 
by a direct gap of about 95 meV. The spectral weight, suppressed by the gap opening, is 
transferred to high energies. It has been long debated up to which energy the total spectral 
weight in FeSi is effectively fully recovered 6 ' 9 . In the data of Degiorgi et al. (Ref. 6), the 
spectral weight is essentially recovered at a frequency uj c ~ 4i? g , therefore without any need 
to invoke an integration of U\ (cu) to very high frequencies. This conclusion is at variance with 
claims 9 , suggesting a redistribution of spectral weight extending up to very high energies, 
even beyond the highest energy limit of that experiment. This feature was also found to be 
indicative for a Kondo insulator description of FeSi (Ref. 9). As far as FeSb 2 is concerned, 
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FIG. 1: FeSb2 sample used for our optical investigations. The crystal axes, determined through 
the Laue spectroscopic experiment, are also shown. 

a comprehensive optical study is still missing. 

In addition to FeSi, FeSb 2 seems to be a promising model system in order to investigate 
the applicability of the Kondo insulator concept with respect to the nearly itinerant magnetic 
semiconductor picture for 3d intermetallic compounds. The absorption spectrum of FeSb2 
should reveal important information about its intrinsic physical properties and should allow 
to extract the phonon spectrum as well as the relevant energy scales, like the hybridization 
gap. To test the applicability of the Kondo insulating concept, it is of interest to investigate 
the redistribution of spectral weight above and below T cr in FeSb2, addressing the issue of 
its conservation. 

We report on our optical investigation on FeSb2- The paper is organized as follows: we 
will first characterize the specimen and describe the experiment, followed by the presentation 
of the optical data. The discussion will first tackle the analysis of the phonon spectrum. 
Particular attention will be then devoted to the issue of the spectral weight distribution. 
The relevance of the Kondo insulator scenario for the description of the metal-insulator 
transition at T cr in FeSb 2 will be emphasized. 
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II. EXPERIMENT AND RESULTS 



Single crystals of FeSb 2 are grown from excess Sb flux, as described in Ref. 1. The crystals 
are cut in a rectangular shape with the long direction corresponding to the b crystalline axis. 
Our sample is shown in Fig. 1, which also displays the crystal orientations. We measured 
the optical reflectivity R{uj) of our FeSb2 specimen, over a broad spectral range (30 to 10 5 
cm -1 ) and as a function of temperature. Light was polarized along the b axis (E || b) and 
along the direction perpendicular (E _L b) to the b axis within the (102) surface (Fig. 1). 
The real part <Ji(uo) of the optical conductivity was obtained through Kramers-Kronig (KK) 
transformation 10 ' 11 of R(u). To this end, we applied standard high-frequency extrapolations 
R(u) ~ lu~ s (with 2 < s < 4) in order to extend the data set above 10 5 cm -1 into the 
electronic continuum. Below 30 cm -1 , R{u) was extrapolated towards zero frequency either 
by using the Hagen-Rubens (HR) law R(oo)=l — 2a/ {uj/adc) for reflectivity data at T > T cr , 
displaying a metallic behavior, or by imposing R{uj < u min )=R(u min = 30 cm -1 ) for data in 
the insulating state at T < T cr . Nevertheless, the overall behaviour of ai{u) is not affected 
by the details of these extrapolations. Details pertaining to the experiment can be found in 
Refs. 10 and 11. 

Figure 2 presents R(cu) and a\{uj) in the far-infrared (FIR) spectral range at 10 and 150 
K for both polarizations of light, while the insets display the spectra at 150 K for both 
quantities over the whole measured spectral range (the reader should pay attention to the 
use of the logarithmic energy scale in the insets). R(u) at 150 K displays a metallic behavior 
for both polarization directions, with a plasma edge onset at ~ 3000 cm -1 (inset of Fig. 2a). 
The anisotropy in R(u) is not very pronounced at 150 K, even though the E _L 6 reflectivity 
is depleted with respect to R(u) along the b axis in the spectral range between 300 and 
1000 cm -1 . With decreasing T and for both polarizations, R{uj) flattens in the far infrared 
range (Fig. 2a), an indication for a metal to insulator (MI) transition. The enhancement of 
phonon-like structures, especially for E _L b, is also observed at 10 K (i.e., in the insulating 
state). R(uj) at 150 and 10 K merge together at ~ 350 and 450 cm -1 , for light polarization 
E || b and E _L 6 , respectively. The anisotropic behavior in the optical response between the 
two polarization directions becomes more evident at low temperatures. In the far infrared 
range, R(u) at 10 K along the b axis is approximately 10% higher with respect to the other 
crystallographic axis. Furthermore, the anisotropic electrodynamic response clearly shows 
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FIG. 2: (color online) a) Optical reflectivity R{uj) and b) real part a\{u) of the optical conductivity 
of FeSbi in the far infrared range, with light polarized along the b crystallographic axis (E \\ b) 
and perpendicular to the b axis (E _L b) at 150 and 10 K. Insets: R(u>) and cr\(uS) for frequencies 
up to 10 5 cm -1 (the reader should note the use of the logarithmic energy scale). Blue and red 
arrows indicate the low temperature phonon modes for E \\ b and E _L b, respectively. 



the presence of different phonon features along the two crystallographic directions. 

The real part (Ji(u) of the optical conductivity, as extracted from the KK transformations, 
is shown in Fig. 2b. At 150 K (inset of Fig. 2b), one can easily recognize the presence 
of a Drude metallic term at low frequencies, for both polarization directions, as well as 
the strong absorption at ~ 20000 cm -1 due to an electronic interband transition. With 
decreasing temperature, the Drude term vanishes along both directions (main panel of Fig. 
2b). Furthermore, for E || b we can clearly distinguish phonon modes at ~ 107, 231, 270 
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cm -1 (blue arrows in Fig. 2b); the latter one seems to be made up by two distinct modes, 
yet not completely resolved. On the other hand, a±(u;) for E _L b displays phonon features 
at ~ 122, 217, 260 cur 1 (red arrows in Fig. 2b). The mode at 217 cm" 1 is clearly present 
at low temperatures only. 

III. DISCUSSION 

A. The fitting procedure 
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The optical conductivity is analyzed within the framework of the phenomenological 
Lorentz-Drude 10 ' 11 model extended to the Fano approach 12 ' 13,14 in order to account for the 
asymmetric shape in the phonon modes. The Drude metallic behavior is defined by the 
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TABLE II: FeSb 2 - E Lb. Lorentz-Drude and Fano fit parameters. All parameters are given in 
units of [cm -1 ], with the exception of eoo and q which is unit-less. 
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formula: 



°d{u) = , (1) 

47T7D 1 - IU/'Jd 



where the plasma frequency is uj p d = a/ AirNe 2 jm and 7d is the scattering rate. Besides 
the Drude term, we add four and two Lorentz harmonic oscillators for E || b and E ± b, 
respectively. The contribution due to a Lorentz harmonic oscillator is: 

vlorentM = ^ ^2 ^ + ~ , (2) 

where uioj is the resonance frequency, jj the damping and u P j the strength of the j-mode. 
For E .Lb, there is additionally a third contribution, described by a Fano lineshape for the 
asymmetric mode at 260 cm -1 . The contribution due to the Fano mode is given by the 
equation 13 ' 14 : 
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Fano 



(u) = ia 0j 




(3) 



with 




(4) 



where the same notation as in eq. (2) has been used. The dimensionless qj Fano parameter 



lineshape is recovered). The Lorentz and Fano contributions describe the modes, marked 
by arrows in Fig. 2b, which clearly develop at low temperatures. 

Table I and II summarize the temperature dependence of the parameters employed in 
our fits. We immediately observe that with decreasing temperature, the plasma frequency 
decreases for both polarizations. From our fits we can define two intervals of temperatures 
for the MI transition, t|^ 7 = 70-100 K and T~mj = 50-70 K, as the range of temperatures 
at which the Drude term vanishes. Our Tmi values are in good agreement with T cr , which 
defines the metal-semiconducting crossover in the p(T) transport data 1 . Another interesting 
aspect in the infrared spectra of FeSb2 is the quite large mode strength of the infrared active 
phonons. Particularly, the phonons at u 02 , ouo3 and c^04 for E || b and the Fano mode at 
cj 3 for E _L b exhibit at low temperatures mode strengths between 600 and 800 cur 1 . This 
could suggest highly ionic character for FeSb 2 - In this respect, FeSb2 shares a common 
property with FeSi (Ref. 9). Schlesinger et al. also pointed out that strong phonons have 
been observed in prototype heavy fermions systems and cuprates, as well 9 . 

While the temperature dependence of the phonon modes and in general of a\(uS) in the 
far infrared (FIR) spectral range will be extensively discussed below, we anticipate here that 
the proposed fit only partially reproduces the measured spectra. Indeed, there is a broad 
continuum of low frequency excitations in the FIR spectral range and at low temperatures 
for both polarizations, which is not encountered by this fit procedure. 

B. Lattice Dynamics 

1. Factor group analysis 



is a measure of the degree of asymmetry of the peak at ujqj (for \q 



oo a Lorentzian 



FeSb 2 crystallizes in the marcasyte-type structure with two formula units per unit cell. 
FeSb 2 is made up of FeSb 6 octahedra which form edge sharing chains along the c axis. The 
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space group is the centrosymmetric Pnnm orthorhombic group 1 , for which we report the 
factor group analysis in Table III. By using the correlation method 15 , and after subtracting 
the acoustic modes, we find for the total irreducible representations: 

Tpnnm = A g + 2A U + B Xg + B lu + 2B 2g + 3B 2u + 2B 3g + 3B 3u . (5) 



TABLE III: Normal modes analysis of FeSb2 for the Pnnm symmetry group. 
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By comparing our data with the theoretical calculation for the number of infrared active 
phonon modes along the b axis direction (E \\ b), one can notice that our optical measurement 
suggests the presence of four instead of three phonon modes. Nevertheless, the resonance 
frequencies of the modes at 257 and 271 cm -1 are quite close to each other. Therefore, 
optical spectra may reveal a slightly lower symmetry for FeSb 2 with respect to the one 
predicted by the Pnnm group 1 . 

Along the direction perpendicular to the b axis (E _L b) three phonon modes are present 
at low temperature, while only two of them are still observable above 50 K. If one assumes 
the Pnnm group symmetry for FeSb 2 , our spectra would imply that the axis perpendicular 
to b in our experiment corresponds to the crystallographic direction a. On the other hand, 
the disappearance of the phonon mode at 216 cm -1 at high temperatures might suggest a 
phase transition into a state, where this specific phonon is silent and not infrared active. 
Nevertheless, no structural changes have been found so far in FeSb 2 down to 10 K (Ref. 2). 
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FIG. 3: (color online) a) Blow-up of <ti(cj) around the mode at 260 cm' 1 at 10 K and for E _L b. 
The fits of this mode with either a Lorentz harmonic oscillator or a Fano lineshape are shown for 
comparison, b) Temperature dependence of the oscillator strength w P 3, c) resonance frequency U03, 
d) scattering rate 73, and e) asymmetry factor q% of the mode at 260 cm -1 within the Fano-like 
fit. The blue dashed line indicates the upper limit of Tj^j. 



2. Fano analysis 

A blow up of the spectral range pertinent to the asymmetric mode at 260 cm -1 for E _L 6 
at 10 K is displayed in Fig. 3a. It is easily recognized that the Fano mode is more appropriate 
than a simple Lorentz oscillator in order to reproduce the overall shape of this mode. The 
temperature dependence of the fitting parameters for this peculiar mode is displayed in Fig. 
3b-e. The oscillator strength is associated to the effective transverse charge and has a non- 
trivial temperature dependence (Fig. 3b) since it first decreases with increasing T, passes 
through a minimum between 50 and 70 K (i.e. in correspondence with the MI transition 
temperature) and then increases again. 

A similar trend is also found for the resonance frequency a;o3 (Fig. 3c), which slowly 
increases with increasing T before displaying a well defined softening at high temperatures. 
X-ray diffraction studies 2 have established that the Fe — Sb — Fe bond angle, associated 
with the edge sharing octahedra along the c axis, increases above 100 K. The T-dependence 
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of the Fe — Sb — Fe angle behaves rather closely to the behavior of uj 3. These observations, 
together with the appearance at low T of the mode at 216 cm -1 , are clear evidence for a 
structural distortion accompanying the MI transition. 

The damping factor 73 (Fig. 3d) increases above T^ fe 7 . Since 7 ~ 1/r , this indicates a 
two times higher phonon lifetime r in the low T insulating state. Such a strong correlation 
between the onset of the insulating state and the loss of phononic relaxational channels also 
suggests the presence of a coupling mechanism between the electronic degrees of freedom 
and the lattice vibrations, as already claimed for FeSi (Ref. 16). 

Finally, Fig. 3e displays the T-dependence of the dimensionless Fano parameter q. At 
high T, \q\ is of the order of 16. Such a high g-value implies that above 100 K an asymmetric 
phonon lineshape can be hardly recognized; a Lorentz harmonic oscillator would fit equally 
well this mode. On the other hand, q takes values between -6.4 and -7.9 below Tj^j, indi- 
cating a sizable asymmetric lineshape. The presence of Fano lineshapes is predicted when a 
discrete lattice vibration is degenerate with states belonging to a continuum which may be 
of electronic or magnetic nature 12 . Interestingly, Fano lineshapes have been also observed 
in the infrared spectra of FeSi (Ref. 16). In that case, the asymmetry was ascribed to 
the interaction of the phonon mode with an electronic resonance peaked in the mid-infrared 
range, i.e. above the phonon resonance frequency. However, for FeSb2 the sign of q(< 0) is 
consistent with the picture of the phonon at ~ 260 cm' 1 coupled to an energy scale which is 
lower than the phonon resonance frequency. The forthcoming analysis of the spectral weight 
distribution in the absorption spectrum of FeSb2 will shed more light on this low energy 
scale. 

C. Spectral weight analysis 

The ai(cu) spectra at 10 K present for both polarizations a deep minimum at about 280 
cm -1 , which coincides with the onset of a broad and temperature independent tail, ascribed 
to the incoherent electronic contribution. We call such an energy, defining the onset of 
the interband transitions, E g . E g ~ 280 cm" 1 is therefore an appropriate estimation for 
the insulating optical gap in FeSb2- The value of E g is sizably larger with respect to its 
estimation based on the activated behaviour of the resistivity, which indicates A p (a, c) ~ 200 
cm -1 , and A p (b) ~ 170 cm -1 (Ref. 1). Similarly to the case of FeSi (Ref. 8), this 
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discrepancy in the gap value among different experiments could be reconciled by assuming a 
small indirect energy gap for the transport properties and a larger direct gap for the optical 
response. On the other hand and at variance with FeSi (Ref. 16), E g does not display 
any temperature dependence up to the temperature (T MI ), at which the MI transition is 
observed. 

By integrating the optical conductivity (i.e., J^ c erf rude {u)duj), we can achieve the spec- 
tral weight encountered in the absorption spectrum over a frequency interval between zero 
and uj c . Because of the temperature dependence of ai(cu) (Fig. 2b), there is an obvious redis- 
tribution of the spectral weight from low to high frequencies in our spectra. A sum rule 10 ' 11 
requires that the integrated area under o\ (uS) must be unchanged at any temperature; there- 
fore predicting the conservation of the spectral weight. In a conventional semiconducting 
scenario, the T-induced MI transition is simply due to the thermal excitation of charge car- 
riers through the energy gap E g . From an optical point of view, one would therefore expect 
that by lowering the temperature through T g =E g /kB the Drude term progressively vanishes 
while the spectral weight lost at uo < E g redistributes just above and close to E g . Such a 
picture is inadequate for FeSbz- We find that the overall depletion of a\{uj) at 10 K leads 
to a removal of spectral weight at low frequencies, which mostly (i.e., 90 %) piles up in the 
spectral range between E g and 3000 cm' 1 . However, a full recovery of the spectral weight 
only occurs when integrating ai(u>) up to 1 eV. The observation, that one needs to go to 
unusually high frequency to satisfy the conduction sum rule, might suggest, similarly to 
FeSi 9 , that the physics of FeSb 2 involves an energy scale much larger than the gap energy 
E g . The transfer of spectral weight in spectroscopies of correlated electron systems has been 
discussed in the framework of the periodic Anderson model for the situation pertinent to 
the Kondo insulators 17 . 

We can now calculate the total spectral weight encountered in cr 1 (a;) of FeSb2 at fre- 
quencies lower than E g . To do so we have subtracted the phonons' contribution. The 
temperature dependence of this latter quantity is plotted in Fig. 4 for both polarization 
directions (open and filled black squares for E \\ b and £16, respectively). We remark 
that at T > Tmi a fraction of the total spectral weight below E g is due to a background of 
weight distributed over an interval of energies extending from up to E g . We associate such 
a background of spectral weight to the presence of a continuum of low lying excitations. 
This continuum of excitations defines a broad far infrared mode (FIM). The total spectral 
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FIG. 4: (color online) Temperature dependence of the spectral weight encountered in g\{uj) at 
u < Eg for E || b (a) and E _L b (b). The black squares represent the integrals of cri(w) up to 
280 cm -1 after having subtracted the fit to the phonon modes. The green triangles represent the 
spectral weight of the metallic Drude term. The blue and red circles represent the spectral weight 
(5A) of the broad infrared mode (FIM) for E \\ b and E _L b, respectively. 

weight up to Eg decreases with decreasing temperature down to T MI (Fig. 4). Below T MI , 
some residual spectral weight (which can not be of vibrational origin) survives inside of the 
gap. As demonstrated in Fig. 4, the spectral weight of the Drude term disappears below 
Tmi- The removed Drude weight is transferred at high energies and partially piles up in the 
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spectral range covered by FIM. The size of this FIM term increases indeed with decreasing 
temperature, as shown in Fig. 4. The continuum of excitations, defining FIM, is better 
highlighted in the insets of Fig. 5. The red and blue area for for E \\ b and E -Lb, respec- 
tively, represent the spectral weight encountered in FIM (for the data at 10 K). Such areas 
(5 A) are calculated by subtracting the spectral weight encountered in the Drude-Lorentz fit 
from the total spectral weight underneath the measured <Ji(u) curve. Figure 5 displays then 
the temperature dependence of 5 A, normalized by its value at 10 K (Ref. 18). At approxi- 
mately T MI the spectral weight associated to FIM suddenly increases and the T-dependence 
5 A vaguely mimics the behaviour of an order parameter. 




FIG. 5: (color online) Temperature dependence of the normalized spectral weight (SA) of FIM. 
Insets show the area 5 A (blue and red for E \\ b and E _L b, respectively) at low temperature, 
encountering the spectral weight of FIM (see text). 



The temperature dependence of 5 A seems to shape the intrinsic electronic properties of 
FeSb2- In fact, the presence of such a (FIM) continuum implies that several localized states 
are present at low temperatures between the valence and the conduction band. Above Tmi 
those states, generating FIM, are thermally occupied. Envisaging an overlapping of these 
states, a coherent electronic transport can then develop. In this context, it is not surprising 
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that Tmi is lower than E g . Furthermore, the presence of this continuum of excitations 
supplies a natural explanation for the Fano behavior of the mode at w 3=260 cm -1 for E Lb. 
As anticipated above, the negative sign of the q Fano parameter implies an interaction of the 
mode at u>03 with a continuum covering an energy spectral range located below uj 03 . Such a 
continuum of excitations is precisely provided by FIM. As comparison, we note that in the 
case of FeSi, where no evidence for an equivalent FIM feature was found, the asymmetry 
of the phonons was on the contrary ascribed to the red-shift of the gap edge 16 . The gap 
edge in FeSb 2 is temperature independent. Therefore, this excludes any possibility for a 
temperature dependent coupling between the phonon modes and mid infrared excitations. 

IV. CONCLUSION 

The electrodynamic response of FeSb2 is characterized by a metal-insulator transition 
occurring at T < 70 K, and opening of a gap E g ~ 280 cm -1 . With decreasing temperature 
there is a depletion of a\{ui) signalling the removal of spectral weight, which is shifted to high 
energies up to 1 eV. Furthermore, we have observed that the suppressed Drude weight with 
decreasing temperature (i.e., at T < T cr ) partially piles up in the broad far infrared mode 
(FIM), which we have ascribed to a continuum of excitations. An ordinary semiconductor 
picture cannot account for the charge dynamic of FeSb2- While some features in the optical 
response may be compatible with a Kondo insulator scenario, it remains to be seen how the 
presence of the continuum of excitations in the far infrared and at u < E g might eventually be 
consistent with such a description. Our findings challenge indeed the theoretical treatment 
of the Kondo insulators within the framework of the periodic Anderson model 17 . 
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